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There have been recent suggestions for searching for the electron electric dipole moment, using 
solid state experiments with compounds containing Gd 3+ ions These experiments could im- 

f*^ | prove the sensitivity compared to present atomic and molecular experiments by several orders of 

^s^j ■ magnitude. The analysis of the problem requires a calculation of the enhancement coefficient K for 

the electron electric dipole moment in the Gd 3+ ion. In this work we perform this calculation. The 
result is K ~ —4.9 ± 1.6. Limitations of the accuracy of the calculation are mainly due to the lack 
. of data on Gd 3+ excitation spectra. We formulate which quantities have to be measured and/or 

calculated to improve the accuracy. 
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I. INTRODUCTION 



A violation of the combined symmetry of charge conjugation (C) and parity (P) has been discovered in the decay 
of the K° meson about 40 years ago ||. The exact origin of this symmetry violation remains an enigma, although 
the standard model of electroweak interactions can describe these processes phenomenologically. It has also been 
proposed by Sakharov Q that the matter-antimatter asymmetry observed in our universe could have arisen from 
a CP-violating interaction active at an early stage of the big bang. The CP-violation implies a time-reversal (T) 
asymmetry and hence violation of the combined TP-symmetry, because there are strong reasons to believe that the 
r"| • combined CPT-symmetry should not be violated ||. An electric dipole moment (EDM) of a system in a stationary 
Oh quantum state indicates a violation of the TP-symmetry; this is why searches for EDMs of elementary particles, 
atoms and molecules are a very important approach to the studies of violations of fundamental symmetries ||. In 
the present work we concentrate on the EDM of the electron d e . 
^ . At present the best limitation on d e comes from the Berkeley experiment with an atomic thallium beam [Q, 
^sO " d e < 1.6 ■ 10 _27 ecm. There are some ideas for improving the sensitivity. One way of improvement is an experiment 
with metastable levels of PbO molecules A breakthrough could be achieved in solid state experiments. This idea 
was suggested by Shapiro in 1968 §). The application of strong electric fields to electrons bound within a solid would 
align the electric dipole moments of these electrons. This should lead to a simultaneous alignment of the electron 
spins; the magnetic field arising from this alignment could be detected experimentally. An experiment of this kind 
has been performed with nickel-zinc ferrite |lO| . However, due to experimental limitations, the result was not very 
impressive. Interest to the approach has been renewed recently due to suggestions of Lamoreaux and Hunter j| to 
perform similar experiments with gadolinium gallium garnet, GdsGasO^, and gadolinium iron garnet, GdaFesO]^, 
employing new experimental techniques. First estimates of sensitivity promise to improve the current upper limit on 
the electron EDM by at least three orders of magnitude, depending on the experimental setup. A thorough analysis 
of the problem requires the calculation of the enhancement coefficient for the electron EDM in the Gd 3+ ion. We 
Oh perform this calculation in the present work. 



II. SINGLE PARTICLE CONTRIBUTION 

The Gd 3+ ion has a nucleus with charge Z — 64 and 61 electrons. The 7 electrons in the outer shell occupy 
the 4/ orbitals. So the shell is half filled, and hence the total orbital angular momentum is zero, L — 0, and total 
spin 5* = 7/2. Let us consider the state with maximum z-projection of the spin. In the representation of second 
quantization the ground state wave function is of the form 

1.9 s ) = /-3T^-2T^-iT^0T^iT^2T^3T 1^)' W 

where f^ na is the creation operator for a 4/ electron with spin a and z-projection of orbital angular momentum l z = m. 
The TP-odd interaction of the electron EDM with the electric field E is of the form, see e.g. Ref. M 
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V d = -d el0 X ■ E, 



(2) 



where 70 and £ = 70757 are Dirac 7-matrices. Because of Scruff's theorem it is crucially important to account for 
very complex many-body screening effects, when working with the Hamiltonian (0); technically this means that the 
many-body perturbation theory practically is not convergent. The standard way |6j to avoid this complication is to 
split the Hamiltonian into two terms: Vd — — d e 7oS • E = — d e T, ■ E — rf e (7o — 1)S • E. Then due to Schiff's theorem 
|pT[ the contribution of the first term to the EDM of the ion is identically zero, so one can reduce the interaction 

V d - VJ = -4(7o - 1)S • E. (3) 

Perturbation theory with this operator is reasonably convergent. 

In leading order of single particle perturbation theory the EDM of the ion D sp is given by diagrams shown in Figj^. 
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FIG. 1. Leading contribution to the Gd 3+ EDM. The dashed line denotes the dipole moment ez — ercosO, and the cross 
denotes the TP- odd reduced interaction VJ , eq (Qj. 



The corresponding formula reads 



D SP = 2£ <4/ mT |er^ g |nKn|^|4/ mT ) = ^ (4) 

mn * n 

where the EDM enhancement coefficient is given by 

v r, \ - ( 4 /mT I {r/a B ) cos 0\n) (n\(ea B / 'd e )V^\Af rn] ) 
mn J 

Here e is the electron charge and as is the Bohr radius. We first consider bd and 5g intermediate states, as the 
contribution of other states is much less. The El matrix elements are of the form 



JQ — m 2 

(4/mTl('"/as)cos6 l |5d mT } = \j — — — r 4fM , (6) 
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(Af m ]\(r/a B ) cos e\5g m ] ) = y — ^ — r 4 /, 5 g, 
where is the El-transition radial integral expressed in atomic units (Bohr radius), 

roc 

r lk - / Ri{r)R k (r)r 3 dr. (7) 
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Ri(r) is the radial wave function of the corresponding orbital, J °° R^(r)r 2 dr = 1. 

The matrix element of VT has been calculated earlier, see e.g. Ref. ||. In jj coupling scheme, and using a 
semiclassical approximation, it reads 

4Z s a 2 Z 2 

(n'jl±l\(ea B /d e M\njl) = ~ ^ _ x) 'Jy^^ ( 8 ) 



where a = 1/137.036 is the fine structure constant, 7 = + 1/2) 2 - (Za) 2 , E = 2Ry = 27.2 cV = 

2 x 109,737 cm -1 is the atomic energy unit, Z e g — 4 is the effective charge for electron motion at large distances, 
and finally the effective principal quantum numbers v and v' are defined by the electron energy with respect to the 
ionization limit, E — — Z 2 s Eq / (2v 2 ) . Our numerical calculations in Hartree-Fock approximation agree with the semi- 
classical result (||) within ten per cent. The ground state wave function (Q) is given in LS coupling scheme, therefore 
we have to rewrite the operator VJ in this representation. A simple calculation gives 
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/q _ m 2 

(5d mT |(ea B /d e )^'|4/ mT ) = -\j <5d 5/2 |(ea B /d e )VT|4/ 5/2 >, (9) 



/ 16 — 777," 

{5g m ^\(ea B /d e )VJ\4f m ^ = -y — — — (5<7 7/2 |(ea B /4)K[|4/ 7/2 ). 

Combining equations (||) , @ , (§) , and @ and performing the summation over m, we find the following expression for 
the single particle contribution to the enhancement coefficient 

K = K7 3 ry 2 7 2 ( V -M^ 1 i 4 r JI^9 1 | nm 

Cff l v 75/ 2 (47 5 2 /2 -l)K/^) 3 / 2 (E 4f -E 5d ) + 3 77/2 (4 77 2 /2 -f)(^ 5ff )3/2 (£7 4/ - ^ 5ff ) ^ " 1 ' 

Note that from here on we skip the atomic energy unit Eq in all equations, assuming that all energies are expressed 
in units of E . To find values of the parameters that appear in this equation, we have performed a Hartree-Fock 
calculation for the Gd 3+ ion. It is known that a Hartree-Fock calculation for an open shell is not a uniquely defined 
procedure. In our calculation we used the following averaging: we assumed that the shell is fully occupied, but the 
occupation number of each single particle /-orbital is 1/2. This is a crude approximation, it gives reasonable values 
of energy levels and radial integrals, but one cannot rely on this calculation as far as energy splittings are concerned. 
As a result of the calculation we obtained the following values of energies and radial integrals: 

E 4f = -1.65 E 5d = -1.20 E 5g = -0.32 (11) 
f4f,5d = 0.63 r 4/i5g = 0.088 

Where available, experimental data should be used, but unfortunately the experimental data on this particular ion is 
scarce. Only the value of the 4/ energy level (ionization limit) is known, see Ref. Jl2| , 

E 4f = -355000cm" 1 = -1.62 E . (12) 



It agrees well with (11). From ( |12[ ) and ( |11| ) one finds v 4 f — 2.22, v^d — 2.58 and v^g = 5.0. 

The most important is the E 4 f — E$d energy splitting. The accuracy of the present Hartree-Fock calculation is not 
sufficient to determine this splitting. There is reliable experimental data for the energy levels of Gd 2+ , Gd 1+ , and 
Eu 2+ , see Ref. Q . Naive extrapolation of the splitting from these ions gives 

A: E 5d - E 4f w 40,000cm" 1 « 0.18 £ - (13) 

On the other hand there is experimental data for Gd 3+ jl3| that indicates 

B : E 5d - E 4f fa 100,000cm" 1 w 0.45 E . (14) 

Unfortunately Ref. (l3) does not contain identifications of all the possible levels, therefore we cannot quite rely on the 
data. For this reason we will present two estimates of K: one for the case (A), see (p"3|), and another for the case (B), 
see (|l4|). Experimental and/or theoretical determination of the E§d — E 4 f splitting would be the most important to 
improve the accuracy of the calculation of the enhancement coefficient K. 

Substituting values of the parameters listed above into eq. (jl^) we find the contribution of the 5d intermediate 
state to the EDM enhancement coefficient: 

A: K 5d = -4.5, (15) 
B: K bd = -1.8 

The contribution of the 5g intermediate state is very small. We also estimate the contribution of the higher d-levels 
(mainly the continuous spectrum) as 

A, B : ^ nd ,„>5«-l. (16) 

Altogether this gives the following value of the single particle contribution to the EDM enhancement coefficient 

A: K sp = -5.5, (17) 
B : K sp = -2.8 
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III. MANY-BODY CORRECTIONS 



In the situation with cesium or any other atom with valent s- or p-electrons || the single particle estimate is 
satisfactory. However here we have /-electrons that have a very small wave function in the vicinity of the nucleus. 
As a result the single particle contribution is strongly suppressed. Technically this suppression is reflected in eq. (||); 
the matrix element is proportional to I/7 3 , and for f 5 / 2 -electrons 7 w 3, while for s 1 / 2 - or p^-electrons 7 « 0.88. 
Therefore it is very important to estimate the many-body corrections. The leading many-body corrections to the EDM 
enhancement coefficient K m \> are shown in diagrams FigJ^, where the wavy line denotes residual Coulomb interaction 

1 °° k Att r k 

vc - - E E ^tt • ^ (is) 

We only account for diagrams with momentum of the Coulomb quantum k not higher than 2. The contribution of 
each diagram from Fig.^J must be doubled, because an opposite order of operators is also possible. 



4f 5pl/2 4f 6pl/2 5s 4f Spl/^/ 
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FIG. 2. Leading many-body corrections to the Gd s+ EDM. The dashed line denotes the dipole moment ercosB, the cross 
denotes the TP-odd reduced interaction VJ, eq (Hj, and the wavy line denotes residual Coulomb interaction 

The formula for the many-body contribution to the enhancement coefficient due to the diagrams in Fig.^ reads 

cos 9\5d m )(5pi/2 1 (5d m | Vc7|6s) |4/ m ) (6s|(eas/ d e )VJ\5pi/ 2 ) 
„ V ( E if ~ E 5d ) (E 5p - E es ) 

| (5s|(ea B /d e )Vy|6pi /2 )(4/ m |(r/a B )cosg|5d m )(6pi /2 K5d m |Vb|5a)|4/ m ) 

(E 5s — Eq p ) (E4f + E 5s — E 5d — E &p ) 
, (4/ m |(r/ aB )cos0|5d„O(5 S |( eaB /4)V^|6p 1/2 )(6p 1/2 |(5d m |Vb|5 S )|4/„ 



(Eif - E bd ) {E if + E 5s - E 5d - E 6p ) 

We have only included s — P1/2 matrix elements of VI and only the intermediate states involving bd electrons which 
give the main contribution. All the diagrams in Fig.g are exchange ones, this is why the sign in eq.(^|) is negative. 
Matrix elements of the Coulomb interactions are of the form 



{5p 1/2 \{5d m \V c \6s)\4f m ) = -i^i-^F( 2 )(5p,4/;5d,6s) (20) 



(6pi /2 |(5d m |Vfc|5s}|4/ m ) = -iyi^!^( 2 )(6p,4/;5d,5 S ) 



where F^ k '(i,j; h, I) is the usual Coulomb radial integral, 



foo j'OC k 

F^(i,j;h,l)= I / -^Rii^Rj^Rh^Ri^rfri^drs, (21) 
Jo Jo r> 

expressed in atomic units. 

Using matrix elements ([?]), (^) and (|2^), and performing the summation over m in (|l^), we obtain the following 
expression for the leading many-body correction to the enhancement coefficient 

8Z 3 q 2 Z 2 ff f ^ 2 )(5p,4/;5d,6 5 ) 1 1__ . . 

mb " 5 7 i/ 2 (4 7l 2 /2 - 1) 4/ ' M I (^ 5P )^ ' (Et, - E 5d ) ■ (E 5p - E es ) [ZZ) 

F( 2 )(6p,4/;5d,5s) 1 /ll 



(v5sV6 P ) 3 / 2 {E±f — E 5d + E 5s — E 6p ) \ (Eif — E 5d ) (E 5s — E, 
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Like the parameters presented in (|LLj) the energy levels and radial integrals in this formula have been calculated using 
a Hartree-Fock method 



E 5s = -3.15 E 5p = -2.26 E 6s = -1.03 E ep = -0.84 
F {2) {bp, 4/; 5d, 6s) = 0.028 F (2) (6p, 4/; 5d, 5s) = -0.023. 



(23) 



The corresponding effective principal quantum numbers are v$ s = 1.59, v§ p — 1.88, vq s — 2.78, v§ v = 3.09. Using 
the E 5 d — E^f energy splitting given in ( ^3[) and (fli}) and substituting all the parameters into eq. (gg), we find the 
many-body correction due to the 6s-, 6p-, and 5d- intermediate states, see Fig. 2 



.4 
B 



K mbl « -0.7, 
ifmbl W -0.3. 



(24) 



One should also perform the summation over higher s-, p-, and d- intermediate states. We estimate this contribution 

as 



mb2 



-0.2 



Combining (24) and (E5J) one finds the leading many-body correction to the enhancement coefficient 



.4 
B 



K mb « -0.9, 
K mh ks -0.5. 



(25) 



(26) 



It is substantially smaller than the single particle contribution ( |l7| ) , so the many-body perturbation theory is conver- 
gent. The final result for the electron EDM enhancement coefficient, K = K sp + K mb , reads 



.4 
B 



K w -6.4, 
K w -3.3. 



(27) 



We recall that the case A corresponds to the energy splitting (13), and the case B corresponds to (J14|). 



IV. CONCLUSION 



We have calculated the electron EDM enhancement coefficient in the Gd 3+ ion. The single particle contribution 
as well as the leading many-body corrections have been taken into account. The result is K w —4.9 ± 1.6. The 
main reason for such a large uncertainty lies in the unknown energy splitting E^f — E^. Experimental and/or 
theoretical determination of the splitting in Gd 3+ would be the most important step for improving the accuracy of 
the enhancement coefficient calculation. 
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